*'  t 


«  3 

'  ►  ■ 


AN  ANALYSIS  OF  THE  TRIAXIAL 
TEST  FOR  COHESIONLESS  SOILS 
by 

Paul  R.  Pas lay  and  Jay  B.  Weidler 


Division  of  Engineering 
Brown  University 
Providence,  Rhode  Island 


May,  1967 


oUW)  1 1 80 

(\i) 


ARPA  SD-86 

Report  E  41,  May  1967 

An  Analysis  of  the  Triaxial  Test  for  Cohesionless  Soils,  P.  R.  Paslay 
and  J.  B.  Weidler 


Errata :  (June  1967) 


1,  Page  13:  Equation  (35)  should  read 

•  „  r  sintfi  3c  ,  , 

°C  =  "2  ^1  -  sin$*3v  +  (c 


.  cos$  3<{u  • 

v 


Page  15:  Equation  (44)  should  read 
S**  =  S**  (v,  o^) 

Z  v.cos2ij(  -  sin<{>  3c  .  \  cos<(> 

=  H*  cos4*  L-2ll  -  sin<|>  3v  U  “  V  (1 gln#)2 

,  tan3ip^C  ~  S^n<^,  , ,cos<J>  [3G  -  2  +  G  cos2if>!k 

sin2<j>*  '  .  2  -  G  [1  +  sin<f>]  ' 

,[1  +  sinj)][3  +  cos2i|Ow3G  2  3G  sin<j>  3c 

2  -  G  [1  +  sin<j>3  ov  3  3a  L1  -  sinij)  3v 


.  {c  -  o.  }  cos<(> 

T  L 


{1  -  sin<{>} 


2  3v 


3.  Page  18:  The  units  of  Equation  (52.c)  should  be 


lb. -in. /lb. 


4.  Page  19:  Equation  (56. b)  should  read 


v  =  -17.48  £ 

n 


5.  Page  19:  Equation  (58.g)  should  read 


I*-  1-755  ! 


6.  Page  32:  Equation  (32)  should  read 


Y  =  -  (tan  40  Z 
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Several  implications  of  an  earlier  analytical  formulation  for 
cohesionless  soils  are  deduced  for  the  triaxial  test  (a  compression  test  in 
the  presence  of  a  hydrostatic  pressure).  A  specific  form  for  the  loading  sur¬ 
face  is  adopted  in  order  to  illustrate  certain  difficulties  associated  with 
uniqueness  and  stability.  In  certain  important  respects  the  limited 
analytical  and  numerical  results  appear  to  agree  with  experience  in  soil 
mechanics  testing. 
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Introduction 

In  a  previous  report  [1]*  the  authors  proposed  an  analytical 
formulation  for  the  prediction  of  isothermal  flow  of  an  idealized  cohesionless 
soil  subjected  to  mechanical  load.  The  material  was  assumed  to  be  time- 
independent,  isotropic,  and  rigid-plastic.  Although,  in  a  broad  sense,  the 
formulation  adequately  represents  certain  aspects  of  the  deformation  character¬ 
istics  of  the  soil,  the  generality  of  the  formulation  must  be  reduced  to  a 
specific  relationship  which  can  be  compared  with  experience  and  experimental 
results . 

Such  a  relationship  not  only  is  useful  in  verifying  the  analytical 
formulation  but  also  provides  a  means  for  evaluating  certain  experimental 
considerations.  In  the  triaxial  test,  for  example,  it  becomes  possible  to 
assess  the  influence  of  the  details  of  the  loading  system,  the  influence  of 
the  membrane  which  encloses  the  sample,  or  the  interplay  between  a  con- 
tinously  distributed  strain  field  and  the  presence  of  slip  planes  on  which 
the  shearing  concentrates. 

The  work  presented  here  is  an  attempt  to  study  some  of  the 
implications  of  the  earlier  analytical  formulation  for  the  specific  case 
of  a  triaxial  test.  This  is  accomplished  by  assuming  a  definite  function 
to  describe  the  loading  surface  and  applying  it  to  a  triaxial  test  on  an 
initially  homogeneous  soil  sample.  This  test  is  currently  the  most  widely 
used  experimental  technique  for  evaluating  soil  response. 

The  next  section  indicates  the  specific  choice  of  the  loading  sur¬ 
face  adopted  in  this  investigation  and  the  stress -deformation  rate 


*  Numbers  in  brackets  designate  references  listed  in  the  bibliography. 
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relations  associated  with  it.  Subsequent  sections  deal  with  questions  of 
stability  and  uniqueness,  and  the  superposition,  onto  a  constant  hydrostatic 
pressure,  of  a  compression  test  which  is  the  triaxial  test.  Finally  calcula¬ 
ted  results  are  given  for  a  specific  case  which  indicate  that  this  formulation 
is  capable  of  predicting  the  type  of  flow  which  is  observed  in  practice. 
Analytical  Formulation 

As  mentioned  earlier  the  material  under  consideration  is  assumed  to 
be  isotropic,  initially  homogeneous  and  to  exhibit  rigid-plastic,  time 
independent,  isothermal  behavior.  The  loading  surface  chosen  is  a  Coulomb- 
Mohr  type  which  closes,  in  some  manner,  at  sufficiently  high  hydrostatic 
pressure.  Figure  1  shows  the  Coulomb-Mohr  surface  in  principal  stress 
space.  A  cross  section  of  this  surface  perpendicular  to  the  vector  i  +  j  + 
k  is  of  interest.  The  plane  of  this  cross  section  is  referred  to  as  the  it 
plane  and  is  shown  in  Figure  2  with  the  notation 


3  <<*!  +  °2  +  °3) 


(1) 


It  is  noted  that  in  the  ensuing  development  tensile  stresses  are  taken  as 
positive.  Figures  1  and  2  are  derived  [2]  by  hypothesizing  that  deformation 
is  initiated  if  the  following  function,  considering  every  plane  through  a 
point. 


resolved  shear  stress |  ^2) 

-normal  stress  +  c 

reaches  a  critical  value  denoted  by  tan  <J>.  In  other  words,  if  considering 
all  planes  through  the  point  in  question,  the  function  (2)  attains  the 
value  tan  <p  on  one  plane  or  simultaneously  or  more  than  one  plane,  i.e. 


resolved  shear  stress 
-normal  stress  +  c 


max 


=  tan  <f> 


(3) 


then  the  associated  stress  state  is  on  the  loading  surface. 
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Clearly,  for  the  Coulomb-Mohr  material,  c  and  <f>  are  the  two  material 

parameters  which  determine  the  loading  surface.  It  has  been  suggested  [3] 

that  such  a  loading  surface  is  only  valid  for  sufficiently  small  values  of 

-a  and  that  for  larger  values  the  surface  closes.  Figure  3  shows  a  section  of 

a  ,  o  ,  a  space  which  includes  the  a.  axis  and  is  parallel  to  the  vector  T  + 

X  2.  o  o 

1  +  k.  The  dotted  line  in  this  figure  is  the  type  of  "cut-off"  of  the  Coulomb- 
Mohr  surface  which  is  envisioned  in  this  work.  The  subsequent  analytical 
development  is  now  restricted  to  be  for  stress  states  on  the  Coulomb-Mohr  por¬ 
tion  of  the  loading  surface.  In  this  way  the  details  of  the  dotted  part  of  the 
yield  surface  in  Figure  3  may  be  ignored. 

In  this  paper  the  symbol  U  will  represent  the  internal  energy  per 
unit  mass  of  the  soil.  It  is  noted  that  in  reference  [1],  U  represented  the 
internal  energy  per  unit  current  volume.  As  in  [1]  U  is  taken  to  be  a  function 
of  the  specific  volume,  v.  Then  following  the  usual  reasoning  that  the  deforma¬ 
tion  rate,  denoted  by  its  principal  values  d^ ,  d d^  and  expressed  in  terms  of 
the  vector  d,T  +  d0]^  +  d  k,  is  determined  or  delimited  by  the  location  of  the 
stress  point  o  ,  o  ,  o  on  the  loading  surface,  restrictions  on  the  direction 

1  4  v 

of  d.i  +  d_3  +  d„k,  may  be  deduced  from  energy  considerations.  Due  to  the  pre- 
X  2.  o 

sence  of  U(v)  the  often  employed  normality  condition  [4]  is  not  directly 
determined.  Instead,  the  acceptable  directions  for  d^T  +  +  dgk,  are  found 

to  be  included  in  a  "fan"  of  directions  at  each  point.  In  order  to  predict 
specific  results  it  is  necessary  to  chose  a  direction  for  d^i  +  d^D  +  d^k.  if 
the  compression  test  superimposed  on  a  hydrostatic  pressure  is  to  be  considered 
a  further  complication  arises  in  that  the  stress  point  is  on  an  edge  of  the 
loading  surface.  The  deformation  rate  chosen  here  which  is  acceptable  according 
to  the  previously  mentioned  derivation  [1]  is  the  sum  of  a  vector  normal  to 
the  loading  surface  and  a  vector  parallel  to  1  +  j  t  k.  As  the  normal  to  the 
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loading  surface  at  an  edge  is  not  uniquely  determined  a  linear  combination  of 
two  vectors  is  determined  for  this  point.  Figure  4  shows  the  details  of  the 
construction.  In  this  case  the  deformation  rate  vector  for  a  triaxial  test  may 
then  be  represented  by 

1  _  f'  l  A 

(4a) 


d,  (A  +  A  )  (  f-.— +  G) 


1_  -  sin_4>_ 
1  +  sin<{> 


d2  =  +  h1  (1  -  G)  +  A2  (  -G) 


d3  =  +  Ax  (  -G)  +  A2  (1  -  G) 


(4b) 

(4c) 


where  A^  and  A^  are  arbitrary  positive  numbers  and 


G  =  a 


9U 

9v 


2  sin<j> 


91J 

3(1  +  sin4>)(a  -  a  +  a  r—  ) 
o  9v 


(5) 


It  should  be  noted  that  when  G  vanishes  the  d^i  +  d^  +  dgk  vector  lies  in  a 
plane  perpendicular  to  the  edge  of  the  loading  surface.  The  advantage  of  the 
formulation  for  G  in  Equation  (5)  is  that  as  a  approaches  oq  for  points  on  the 
loading  surface  the  deformation  rate  vector  approaches  being  parallel  to  the  ir 
plane  (i.e.  -  approaches  constant  volume  flow).  This  is  consistent  with  ex¬ 
perience  in  soil  mechanics. 

Finally  for  a  soil,  similar  to  a  metal,  the  yield  surface  changes  as 
the  deformation  ensues.  The  most  important  variable  controlling  the  changes 
appears  to  be  the  specific  volume  so  that  in  the  following  it  is  assumed  that 
the  material  parameters  c,  <j>,  aQ,  and  U  are  each  function  of  v.  Restrictions 
on  the  behavior  of  these  parameters  are  given  in  Equations  (6)  below.  Although 
the  physical  significance  of  these  parameters  has  been  indicated  in  Figures  1 
through  4  a  brief  discussion  of  them  is  included  for  clarity.  The  parameter 
c(v)  is  a  measure  of  the  cohesive  strength  of  the  material  which  for  truly 
cohesionless  soils  (i.e.  dry  sand  grains)  may  be  assumed  equal  to  zero.  <J>(v) 
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is  often  called  the  angle  of  internal  friction  of  the  material  and  in  this  work 
serves  as  a  measure  of  the  linear  rate  at  which  the  Coulomb-Mohr  surface  diverges 
from  the  hydrostatic  axis  with  increasing  hydrostatic  pressure.  The  function 
ao  (v)  represents  the  spherical  portion  of  a  stress  tensor  on  the  loading  sur¬ 
face  which  would  induce  constant  volume  flow.  Finally  the  parameter  U(v)  is 
the  internal  energy  per  unit  mass. 

The  following  restrictions  on  the  behavior  of  these  parameters  appear 
to  be  physically  appropriate 


c(v)  >  0 

(6a) 

p-  <  0 

9v  — 

(6b) 

0  <_  <(>(v)  <  ir/2 

(6c) 

f£  <  0 

9v  — 

(6d) 

0  <  0 

0  9v  — 

(6e) 

9a 

(6f ) 

1  £  1  <  «  1  -O  1 

(6g) 

The  last  of  these  restrictions  coupled  with  restriction  (6e)  indicates  that 

even  though  internal  energy  may  be  returned  to  the  system  during  plastic 

work  with  increasing  specific  volume ,  its  rate  of  return  with  respect  to  the 

3U 

increasing  specific  volume  is  small  compared  to  a  .  That  portion  of  which 

is  available  to  assist  the  superimposed  stresses  in  doing  plastic  work  while 

3U 

inducing  a  differential  change  of  plastic  strain  is  denoted  by  a  -5—  where  a, 

oV 

a  material  parameter,  is  restricted  to  0  <  o  <  1. 

Uniqueness  and  Stability 

It  has  been  previously  shown  [5]  that  normality  of  the  deformation 
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rate  vector  to  the  yield  surface  must  exist  to  ensure  certain  uniqueness 
properties.  Since  normality  is  not  satisfied  here  a  general  lack  of  uniqueness 
is  to  be  expected.  The  question  arises  then  as  to  how  the  actual  solution 
can  be  selected  from  all  those  deduced  by  the  theory.  That  such  a  selection 
should  be  possible  is  a  result  of  the  observation  that  experiments  in  soil 
mechanics  are  repeatable.  What  is  desired,  then,  is  a  means  of  selecting  from 
the  predicted  solutions  that  one  which  will  occur  in  an  experiment.  In  this 
investigation  of  the  triaxial  test  a  lack  of  uniqueness  of  the  flow  field  is 
shown  and  the  actual  field  must  be  chosen. 

The  criterion  which  will  be  used  to  discriminate  between  solutions 
is  determined  from  the  way  in  which  the  rate  of  work  changes  during  a  pre¬ 
scribed  flow.  Suppose  that  the  external  loading  necessary  to  initiate  the 
flow  has  a  distribution  represented  by  F  and  that  the  associated  velocity 
field  is  represented  by  X.  Now  for  certain  parts  of  the  flow  field  (e.g. 
possibly  part  of  the  boundary)  the  X  vector  is  prescribed  as  a  function  of 

time.  After  the  first  interval  of  loading  during  a  time  interval  At  the 

-> 

displacements  are,  to  first  order  terms  in  At,  XAt.  During  this  period  the 
loading  distribution  is  changed  to  F  +  FAt.  Consider  the  quantity 

■> 

F  .  X  (7) 

system 

which  is  the  time  rate  of  change  of  doing  work  on  the  system.  For  each 
solution  deduced  from  the  theory.  Expression  (7)  can  be  evaluated.  It  is 
hypothesized  here  that  of  all  the  possible  solutions  the  one  which  minimizes 
Expression  (7)  is  the  actual  one. 
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Triaxial  Test  of  an  Initially  Homogeneous  Sample 

Figure  5  shows  a  sketch  of  the  usual  triaxial  test,  a  represents 

L 

the  hydrostatic  pressure  which  is  assumed  to  be  applied  initially  and  main¬ 
tained  at  a  constant  predetermined  value.  For  soils  o  is  usually  a  negative 
number.  Subsequently  a  compressive  stress,  -cr^,  is  applied  to  the  homogeneous 
sample  until  deformation  occurs.  In  this  case  take 

°1  =  °C  +  °L  (8a) 

a3  =  a2  =  oL  (8b) 

From  Figure  2  these  stresses  are  substituted  into  the  stress  state  corresponding 
to  point  A  with  the  result  that  the  value  of  a required  to  initiate  flow  is 


o  =  -  2(c  -  a.)  j— — ■ (9) 

C  LI-  sm<t> 

The  deformation  rates  associated  with  this  stress  state  are  now  investigated. 
Two  separate  solutions  will  be  presented  and  examined.  Alternatively  a  linear 
combination  of  the  two  solutions  could  be  investigated  with  the  same  result. 


In  any  event  the  deformation  rates  are  deduced  from  Equations  (4)  to  be 
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A,  =  A2  >  0 

and  Equation  (10.2)  becomes 

d  =  -  2A_  (  y  4-n4  +  G) 
z  11  +  sm<j> 

Now  rather  than  use  it  is  more  convenient  to  work  in  terms  of  the 

.  _  2 

stant  Z.  Since  d  equals,  in  this  case,  =-,  Equation  (12)  yields 

Z  Li 


.  _  1  ,  Z  ,  ,  _ 1_  +  sin<{> _ 

1  2  Z  1  -  sin4>  +  G  [1  +  sin<j>] 


) 


Now  applying  Equations  (11)  and  (13)  to  Equations  (10)  gives 


,  _  ,  _  _  1  r  (1  +  sin<fi )  (1  -  2G)  ,,  Z, 

x  y  2  1  -  sin<Ji  +  G  (1  +  sin4>)  Z 

In  order  to  evaluate  Equation  (7)  note  that 


F.X=(a  +o''  itR2  Z  +  a  2itRZR 
t  l  L» 

system 


Recalling  that  Z  and  a  are  constant.  Equation  (15)  yields 

Li 


a  < 

dt 


F  .  X  )  =  a  iR^Z  +  (o  +  o  )  2irRRZ  +  o  2itZ(R2  +  RR  ) 
C  CL  L 

system 


The  evaluation  of  a  is  accomplished  by  differentiating  a  in  Equatic 

•  ,  »  ( c  -  oT  )  cos<j>  -  , 

o  =  -2  [  (  -S-^ _  )  *£■  +  - _ - _ —  ]  V 

C  1  v  1  -  sin<£  ;  8v  .  a,2  9v  j 

(1  -  sm<{>) 


Now  v  is  determined  as  follows 

*  <a  ^  j  ^  j  \  r2  sin$  -  3G  (1  +  sin4>)  n  /  Z  , 

v  =  (d  +  d  +  d  )  v  =  -  [- - nTT-rr: — -r-’-rvJ  (  7  )  v 

x  y  z  1  -  sm$  +  G  (1  +  simf)  Z 


(ID 

(12) 

con- 

(13) 

(14a) 

(14b) 

(15) 

(16) 

n  (9) 
(17) 


(18) 
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The  quantity  R  is  determined  by  noting  that  d^  for  this  case  is  equal  to  (R/R) 
which  from  Equation  (14. b)  becomes 

n  _  1  r  (1  +  sin<j>)(l  -  2G)  t  „  Z 

R  '  "  2  L  (l"''Yin*)  J  R  *  Z  U9; 

Differentiation  with  respect  to  time  of  Equation  (19)  gives  for  R 

*  ■  -  -i R  <i>  “r^rw  ” a  frSW*  (<i>  -  <§» 

2  (1  -  2G)  cosiji  |^-  -  (1  +  sin<(>)(3  -  sin^i)  ~- 
'  i  9v  dv, 

+  v  { - 5 - ' 

(1  -  sin<{>  +  G  [1  +  sin$]) 

a  (1  +  simj»)(3  -  sin<)>)  t— 

-  ( - — )  ]  (20) 

"  (1  -  sin<j>  +  G  [1  t  sin4> 3 ) 

•  • 

Now  Equations  (9),  (17),  (19)  and  (20)  may  be  used  to  eliminate  R,  and 

R  from  Equation  (16).  Subsequent  to  the  operation  eliminate  v  with  Equation 
(18).  The  result  may  be  expressed  as 


F  .  X  )  =  ttR2Z  (|)2  S* 


system 


where 


S*  =  S*  (a  ,  v) 
Li 


-  +  2  [{(c  ~  V, !!!!  _ 

(1  -  sin*)2  9V  1  -  sin* 


tl  sin*  -  3G  (1  t  sin*)  }  ]  v  +  2(c  -  o  )  [{s^~ _ } 

*(1  -  sin*)  +  G  (1  +  s in )  L  1  -  sin* 


,(1  +  sinjMl  -  2G)  n  r_  f(l  +  sin<j»)(l  -  2G)  , 

*1  -  sin*  +  G  (1  +  sin*;  °L  ^1  -  sin*  +  G  (1  +  sin*) 
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1  r(l  +  sinjoa  -  2G)  ,2  f(l  +  sinfKl  -  2G)  . 

2  1  -  sin4>  +  G  (1  +  sin<J>)  x  -  sinif>  +  G  (1  +  sin<j>) 

(1  L  (1  +  sin(j))(l  -  2G)  ,  ,  2  sin<fr  -  3G  (1  +  sinjih 

2  1-  sin<}>  +  G  (1  +  sin4>)  1  -  sin^  +  G  (1  +  sin<£ 


2  (1  -  2G)  cosif  -  (1  +  sin<{>)(3  -  sin<f>) 

r  o  V  oV 


3G 


(1  -  sin<£  +  G  [1  +  sin<J>])‘ 


.  (c  -  or  )  COS $  .  , 

2  t  L  3 4>  s  in  $  u\^  ■, 

3  ~~(1  -  sin4>)^  9v  r  -  sin4>  9v 


9c 


.  v  .  {- 


(1  +  sin<f>)(3  -  sinij)) 


9G 
3  a 


(1  -  sin<j>  +  G  [1  +  sin<Ji])‘ 


}] 


,2  sinij)  -  3G  (1  +  sintft)  ■, 
*1  -  sin<}>  +  G  (1  +  sini)>r 


(22) 


S*  may  be  regarded  as  a  function,  for  a  specific  material,  of  v  and  a  since  <p, 
c,  U,  oq  and  consequently  through  Equation  (5),  G  are  considered  to  be  specified 
functions  of  v.  When  S*  (a  ,  v)  is  negative  then  the  system  is  unstable  in  the 
sense  of  [5]. 

For  the  second  deformation  solution  of  the  triaxial  test  a  combination 
of  shear  and  lateral  expansion  of  an  inclined  layer  of  the  material  is  consider¬ 
ed.  The  thickness  of  the  layer  is  taken  as  H.  Figure  6a  indicates  a  plane 

parallel  to  the  layer  and  Figure  6b  indicates  the  deformation  state  after  a 

•  •  *  •  • 

small  amount  of  motion.  The  velocity  field  with  components  u^ ,  u^,  u^  in  the 

deforming  region  may  be  expressed  by 


Un  =  H  ? 


(23. a) 

(23. b) 


UC  =  ° 


(23. c) 


and  are  velocity  components  of  the  upper  rigid  part  of  the  sample  shown 
in  Figure  6b  which  undergoes  a  rigid  body  translation.  The  velocity  components 


TEOUTH- H  at,  I.T3L ARY 

•p --  T\  'i_  F  ’ 

. 1 

"ill 
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Z  and  Y  are  related  to  V  and  V„  as  follows 

n  5 


Z  =  V_  cosiJ>  -  V  sinii 

£  n 


(24. a) 


Y  =  V_  sinil/  +  V  cosifi 

£  n 


(24. b) 


For  this  motion  the  components  of  the  deformation  rate  tensor  in  the  deforming 
region  are 


=  H  V£ 


(25. a) 


d  =0 

nn 


(25. b) 


d  =0 
CC 


(25. c) 


d5n  =  2H  VC 


(25. d) 


d  =  0 

n£ 


(25. e) 


=  ° 


(25. f) 


To  determine  how  these  deformation  rate  components  may  be  matched  with  those 
given  by  Equations  (10)  it  is  noted  that  the  c  and  x  axes  are  coincident,  thus 
d^  must  also  vanish  and  this  requirement  leads  to 

a2  =  r~G  Ai  (26) 

The  principal  deformation  rate  components  may  then  be  expressed  as 


d  =  0 
x 

dy  -  Ai  (rTr} 

d  =  -  A  (  r^-p)(rr^T  +  G) 
z  1  1  -  G  1  +  sin| 


(27. a) 


(27. b) 


(27. c) 


Now  the  angle  <p  in  Figure  6  is  determined  from  the  condition  that  d  be  zero. 

The  Mohr  strain  circle  for  this  case  is  shown  in  Figure  7.  In  order  for  d  to 

nn 

vanish  it  may  be  deduced  from  the  circle  that 
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cos  2<J/ 


and  also  that 


-2  sinij>  +  3G  (1  +  sin$) 
2  -  G  Cl  +  sXn<j>) 


=  Ai  1 


j-2  sin<j>  -  3G  (1  +  simfi)-. 

L  _____  j 


d  =  A,  1  [-1--  sln2* 

n?  1  1  -  G  1  +  sm| 

Combining  this  result  with  Equations  (25)  gives 

H 


VA1  1 


2H 


,-2  sirup  -  3G  (1  +  sin<j>)-. 
1  +  sin4> 


fl  -  (G/2)(l  +  sintj) ) .. 

L"  T  •  — .  J _ X  J 


n  1  1  -  G  u  1  +  sin<|) 
and  introducing  Equations  (30)  into  Equations  (24)  yields 


Z  =  A 


1  1  -  G 


H  rf2  sin<j)  -  3G  (1  +  sind>)->  , 

[{ - T  -t  ■-•"-t - COSlf) 

1  +  sinq 


f2  -  G  (1  +  sim}>),  .  a  .  •  .  -i 

"  { - TTil - }  sin  2 *  sln*] 


(28) 

(29. a) 

(29. b) 


(30. a) 


(30. b) 


(31. a) 


*  .  H  r/2  sin<f>  -  3G  (1  +  sin<j>)i  .  , 

Y  *  h  i—  [< - nrns - -> sin* 

+  {- - y  1"  t-  sin  cosijO  (31. b) 

l  +  sin 9 

Expressing  Y  as  a  function  of  Z  by  use  of  Equations  (31)  gives  after  some 
manipulat ion 

Y  =  (tan  3ip)Z  (32) 

Finally  the  rate  of  change  of  specific  volume,  defined  by  (d^  +  d^  +  dz)  v,  is 


v.  cos  2i l>  Z 


(33) 


cosi|/  '  H 

The  stress  state  change  is  more  complex  for  this  motion  than  it  was  for  the 
first  deformation  mode  considered.  In  the  present  case  the  principal  stress 
directions  rotate  once  motion  begins  and  this  must  be  accounted  for  in  determining 
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-ST*  (  F  .  X).  That  the  principal  stress  directions  must  rotate  is 

J  system 

explained  by  first  noting  that  the  existence  of  v  implies  the  existence  of  $ 

e 

and  G.  Since  Equation  (28)  shows  that  ij)  =  tfi  (<pt  G),  then  if  <f>  and  G  change  dur¬ 
ing  the  motion,  i p  may  also  change.  But  if  the  same  material  which  initially 
deforms  at  the  onset  of  motion  is  to  be  the  only  material  deforming  during  subse- 
quent  motion,  then  =  0.  The  only  way  this  can  occur  is  for  the  principal 
stress  directions  to  rotate.  Associated  with  this  rotation  rate  is  the  develop¬ 
ment  of  a  shear  stress  rate,  o^.  During  an  actual  test  this  shear  stress  rate 
would  have  to  be  sustained  by  the  testing  apparatus.  In  the  following  the 
apparatus  is  assumed  capable  of  adopting  to  the  required  loading. 

For  purposes  of  illustration  consider  the  change  in  the  principal 
stress  directions  that  occurs  in  a  small  time  interval  At  after  motion  begins. 
Retaining  only  first  order  terms  in  At,  the  Mohr  stress  circle  for  the  yz  plane 
is  shown  in  Figure  8.  The  principal  stresses  after  an  increment  of  time.  At, 
are  consequently  °c  +  °L  +  °L*  °L*  rate  at  which  the  principal 

stresses  rotate  is 

a  =  (34) 

CTC 

The  condition  that  these  principal  stresses  continue  to  satisfy  the  yield 
criterion  is  the  same  as  given  by  Equation  (17)  or 


a  =  -2  —  +  (c  -  a  ) 

C  L1  -  cos  $  3v  v  1/ 


cos4>  3$-,  • 

[1  -  sin<f>)^  3v  V 


By  differentiating  Equation  (28)  with  respect  to  time  and  subsequently  setting 


equal  to  -  a  from  Equation  (34)  leads  to 

-&■  si„2*  =  +  .  <3G  -  2.  G 'COS?* ) 

<?£  2  -  G  (1  +  sinifi)  1 

,  a  +  sin<f>)(3  +  cos2t l>)  n 
+  G  [ — 2---GTi“T'VinT5'  ] 


(36) 
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The  terms  <f>,  U,  and  oq  are  functions  of  v  so  that  from  Equation  (5)  it  is  seen 

that  G  may  be  regarded  as  a  function  of  v  and  a  or  for  the  triaxial  test  v 
°C 

and  a  +  —  so  that 

L  u 


•  3G  *1  3G 
G  =  Fv  v  +  3  3?  °C 


where 


3G  32U 
3v  =  a ;  2  • 


2  sin<J>  (a  -  oj 3  -  a  ) 

O  C  L 


3v^  3  (1  +  sin<j>)(cr  -  a  / 3  -  a  +  a  -|^-)2 

O  C  Ij  dV 


(37) 


3a 

c 

3v 


3U  _ _ 

3y  3U,  2 

3(1  +  sin<f>)(o  -  a  / 3  -  a  +  a  — ) 
o  L  L  dV 


3*  3U 

+  — a  — 
3v  3v 


2  cos<)> 


0 

3(1  +  sin$)  (a  -  a  /3  -  a  +  a  r— ) 
o  C  L  3v 


(38a) 


3G  3U  0  . 

r—  =  a  —  .  2  sm<}> 

3v  3v  — - z - - - gp—A 

3(1  +  sin<j>)(a  -  a  /3  -  aT  +  r— ) 

o  C  L  3a 


(38b) 


and 


*3 


3v 


.  v 


(39) 


The  external  loading  in  this  case  does  work  on  the  system  at  a  rate  given  by 


-*  -t  2  2  L  v 

F  .  X  =  irR  o  Z  +  irR  H  -=~r  - 
C  cos^  v 


(40) 


system 


The  quantities  a  ,  H,  ip,  and  Z  in  Equation  (40)  are  constant  and  hence  the 
L 


measure  which  determines  the  actual  flow  field  becomes 


2-  •  2  °L  v  2* 

F  .  X  =  irR  a_Z  +  ttR  H  -=—  (  -  )  +  ttR  a  Y  (41) 

C  cosip  v  yz 

system 

The  third  term  is  present  on  the  right  hand  side  of  Equation  (41)  since  the 
principal  stress  axes  rotate  once  motion  begins.  All  of  the  quantities  appear- 


d_ 

dt 
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ing  in  Equation  (41)  have  been  evaluated  except  (v/v)  which  may  be  found  by 
differentiating  Equation  (33)  divided  by  v  with  respect  to  time  to  obtain 

(  ~  )  =  0  (42) 

Equations  (35),  (42),  (36)  and  (32)  may  now  be  substituted  into  Equation  (41) 
to  obtain 
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are  shown  in  Figure  9.  It  should  be  noted  also  that  when  S  is  negative  and 
H  approaches  zero,  S  approaches  an  unbounded  negative  number.  As  no  parameter 


exists  in  the  first  solution  such  as  H,  which  is  arbitrary  (0  <_  H  <_  Z  cos*J>  - 
*  ** 

2R  sinijO ,  S  is  bounded.  Consequently  when  S  is  negative  the  solution  depicted 

**  ... 

in  Figure  9b  must  be  the  actual  solution.  When  S  is  positive  then  its  minimum 


Value  must  be  compared  to  S  to  determine  the  actual  mode  of  deformation. 


Clearly,  the  above  solutions  may  not  be  the  only  solutions  so  that  the 


comparisons  given  here  really  only  pertain  to  determining  which  which  of  the 


considered  solutions  is  likely  to  be  the  actual  solution. 


Numerical  Solution 

In  order  to  demonstrate  predictions  which  the  preceding  results  yield 
a  numerical  example  is  presented.  Several  functions  must  be  given  explicitly  in 
order  to  make  this  evaluation.  Published  experimental  results  [6]  for  Leighton 
Buzzard  sand  were  used  in  order  to  select  appropriate  input  for  the  calculations. 
The  following  choices  are  consistent  with  the  experimental  results  and  theoretical 
requirements  of  Equations  (6): 

c  =  0  (47. a) 


*  =  *o  +  to 


a 

o 


fv 


(  v  -  v  ) 
o 

V 


(47. b) 
(47. c) 


aU  =  aK 


a  dv 


(47. d) 


where 


34> 

3v 


>  0 

<  0 


(48. a) 


(48. b) 
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V  <  V  <  V, 

o  1 

(48. c) 

M  >  0 

(48. d) 

1  >  K  >  0 

(48. e) 

0  <  a  <  1 

(48. f) 

From  Equations  (47)  and  (48)  it  can  be  determined  that 

v,  -  v 

aU  =  aKM  f  -  (v  -  v)  +  (v  -  v  )Jn  ( - -)  ] 

o 

(49. a) 

9U  „„  /l  ‘ 

°  17  *  -  “KH  (v  -  v  > 

o 

(49. b) 

-2  (v  -  v  ) 

a  ^4  =  «KH  1  ° 

9v^  (v  -  v  y 

o 

(49. c) 

9a  (v  -  v  ) 

°  _  .  1  o 

3v  +  /  s2  M 

(v  -  v  ) 
o 

(49. d) 

From  the  experimental  results  of  [6]  the  following  magnitudes  were  selected  for 
the  material  parameters  of  Equations  (48) 


<j>  =  0.686  radians  =  39.4° 

o 

(50. a) 

-  0.100  radians  -  lb. /in.3 

(50. b) 

v  =  16.08  in.3/lb. 
o 

(50. c) 

V;L  =  19.00  in.3/lb. 

(50. d) 

M  =  129  lbs. /in. 2 

(50. e) 

aK  =  0.1 

(50. f) 

The  initial  homogeneous  state  of  the  sample 

for  the  calculations  of  the  numerical 

example  is  described  by 

v  =  16.88  in.3/lb. 

(51. a) 

o.  =  -  50  lbs. /in. 2 

L 

(51. b) 
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The  preliminary  evaluations  needed  for  determination  of  both  S*  and  are 


<f>  =  34.7° 

from  Eq.  (47.b) 

(52. a) 

o  =  -343  ]bs./in.2 

0 

from  Eq.  (47.c) 

(52. b) 

aU  =  +  21.5  lb. -in. /in. 

from  Eq.  (49. a) 

(52. c) 

2  5 

=  +  592  lb.  /in. 

from  Eq.  (49.d) 

(52. d) 

a  ~  =  -  34.3  lb. /in. 2 

from  Eq.  (49.b) 

(52. e) 

2 

3  II  9  r 

o  =  59.2  lb.  /in. 

from  Eq.  (49.c) 

(52. f) 

oc  =  -132  lbs. /in. 2 

from  Eqs.  (9),  (47. a) 

(52. a) 

(52. g) 

a  =  -94  lbs. /in. 2 

from  Eqs.  (1),  (52. g) 

(52. h) 

G  =  +  0.0293 

from  Eqs.  (5),  (52.e) 

(52. b),  (52. h) 

(52. i) 

=  +  0.0140  lb. /in. 3 

from  Eqs.  (38. a),  (52. f) 
(52. b),  (52. g),  (52. e) 
(52. d) 

(52. j) 

■gjj-  =  -  1.038  x  10“  in.  /lb. 

from  Eqs.  (38. b),  (52. c) 
(52. b),  (52. g) 

(52. k) 

For  the  first  mode  of  deformation  considered 

• 

v  =  -  35.4  (|) 

from  Eq.  (18) 

(53. a) 

os-  1570  (§.) 
c  z 

from  Eq.  (17) 

(53. b) 

• 

R  =  -  1.55  R  (|) 

from  Eq.  (19) 

(53. c) 

••  *7  o 

R  =  -  11.78  R  (|0 

from  Eq.  (20) 

(53. d) 

Substituting  the  above  into  Equation 

(16)  yields 
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d_ 

dt 

r 

?  .  X  .  ,r2z 

(Z/Z)2  (+  88) 

(54) 

system 

that  from  Equation  (21) 

* 

S  = 

:  +  88 

(55) 

For  the  second  mode 

of  deformation  considered: 

4>  = 

60.4° 

from  Eq .  (28) 

(56. a) 

• 

v  = 

-  17.48  | 

from  Eq.  (33) 

(56. b) 

. 

a  = 
c 

• 

-  774  1 

from  Eq.  (35) 

(56. c) 

. 

♦  = 

• 

+  1.748  | 

from  Eq.  (39) 

(56. d) 

* 

G  = 

• 

-  °-217  1 

from  Eq.  (37) 

(56. e) 

• 

a 

yz 

from  Eq.  (36) 

(57. f) 

• 

Y  = 

-  1.755  | 

from  Eq.  (32) 

(58. g) 

Substituting  Equations  (42)  and  (56)  into  Equation  (41)  yields 


F  .  X  =  itR2Z  (Z/Z)2  (  -  1023  §  )  (59) 

n 

system 

so  that  from  Equation  (43) 

S**  =  -  1023  (§•)  (60) 

n 

From  the  remarks  previously  made  concerning  the  effect  of  the  sign  of 
S  ,  the  actual  mode  of  deformation  of  the  specimen  in  this  numerical  example  is 
the  one  shown  in  Figure  9b.  There  is  of  course  a  geometric  requirement  that  the 
height  Z  of  the  sample  be  large  enough  to  allow  a  slip  plane  to  form.  That  is, 
for  the  second  mode  of  deformation  to  occur 


z 

—  >  2  tantj/  =  3.52 

K 
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(61) 

It  may  be  noted  from  the  calculations  that  the  first  mode  of  deformation  is 

is 

stable  as  S  is  positive,  primarily  because  the  rate  at  which  the  radius  in¬ 
creases  is  enough  to  offset  the  decreasing  a  . 

c 

The  above  numerical  results  indicate  a  behavior  which  is  consistent 
with  observed  soil  tests.  The  reader  should  note  that  while  this  theory  pre¬ 
dicts  a  flow  similar  to  that  shown  in  Figure  9b,  it  must  be  noted  that,  in 
practice,  a  rubber  membrane  surrounds  the  sample  and  that  the  finite  size  of 
the  grains  prevents  H  from  approaching  zero.  Such  stabilizing  effects  may  be 
responsible  for  the  observed  experimental  result  that  the  flow  shifts  from 
plane  to  plane  as  the  motion  proceeds. 
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Coulomb-Mohr  Loading  Surface  in  Principal  Stress  Space 
Intersection  of  Coulomb-Mohr  Loading  Surface  with  it  Plane 
for  a  =  j  (c  +  +  a3>  <_  c 

Intersection  of  Loading  Surface  and  Plane  in  Principal  Stress 
Space  Which  Includes  the  axis  and  is  Parallel  to  the  Vector 
I  +  j  +  k 

Same  Section  as  Figure  3  showing  direction  for  Deformation 

-►  «► 

Rate  Vector,  d^i  +  d2j  +  d^k 

■>  •> 

it  Plane  Showing  the  Fan  of  Possible  Projections  of  d^i  +  d2j 
+  dgk  onto  the  Plane 

Triaxial  Test  Sample  Showing  Applied  Stresses  o  a  and 
Dimensions 

Deformation  of  Shearing  Layer 

Mohr  Strain  Circle  for  Shearing  Layer 

Mohr  Stress  Circle  for  Time  t  +  At  in  the  yz  Plane  for  Shearing 
Layer 

Modes  of  Deformation  for  Shearing  Solution 


FIG.  I  COULOMB- MOHR  LOADING  SURFACE  IN 
PRINCIPAL  STRESS  SPACE 


FIG. 2  INTERSECTION  OF  COULOMB-MOHR  LOADING 

SURFACE  WITH  7 T  PLANE  FOR 
I 


FIG. 3  INTERSECTION  OF  LOADING  SURFACE  AND 

PLANE  IN  PRINCIPAL  STRESS  SPACE  WHICH 
INCLUDES  THE  0"3  AXIS  AND  IS  PARALLEL 
TO  THE  VECTOR  T  +  J+~k 


(b) 

FIG.  4  (a)  SAME  AS  FIG.  3  SHOWING  DIRECTION  FOR 


DEFORMATION  RATE  VECTOR  ,  i  d,  +  jd2  +  k  d3 , 
WHEN  $  ,<T0  <0- 

(b)  -77*  PLANE  SHOWING  THE  FAN  OF  POSSIBLE 
PROJECTIONS  OF  Td,  +~Jdz  +~kd3  ONTO  THE 
PLANE. 


CTC  +  CTL=  AXIAL  STRESS 


Z 


\/  \ 


RIGHT  CIRCULAR  CYLINDER 
RADIUS  =  R 


CTl  =  LATERAL  STRESS 

TENSILE  STRESS  >0 

crL<  0 

c'c  <  o 


FIG. 5  TRIAXIAL  TEST  SAMPLE  SHOWING  APPLIED 
STRESSES  CTC  ,  crL  AND  DIMENSIONS. 


FIG.  6  DEFORMATION  OF  SHEARING  LAYER. 


SHEAR 


d^£) 

FIG. 7  MOHR  STRAIN  CIRCLE  FOR  SHEARING  LAYER 


SHEAR 


FIG. 8  MOHR  STRESS  CIRCLE  FOR  TIME  t+At  IN 
THE  Af,f  PLANE  FOR  SHEARING  LAYER. 


s**>0 


(a) 


s  **<  0 


(b) 


FIG.  9  MODES  OF  DEFORMATION  FOR  SHEARING 
SOLUTION. 


